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We investigate the folding and forced-unbinding transitions of adsorbed semiflex- 
ible polymer chains using theory and simulations. These processes describe, at an 



elementary level, a number of biologically relevant phenomena that include adhesive 



interactions between proteins and tethering of receptors to cell walls. The binding 
interface is modeled as a solid surface and the worm-like chain (WLC) is used for 



the semiflexible chain (SC). Using Langevin simulations, in the overdamped limit, we 

in 

On 
O 

interactions, temperature, and the persistence length l p we obtained the monomer 



examine the ordering kinetics of racquet-like and toroidal structures in the presence 
of attractive interaction between the surface and the polymer chain. For a range of 



density distribution n(x) (x is the perpendicular distance of a tagged chain end from 
the surface) for all the relevant morphologies. There is a single peak in n(x) inside 
the range of attractive forces b for chains in the extended conformations while in 

o" 

O ■ racquet and toroidal structures there is an additional peak at x ~ b. The simu- 



lated results for n{x) are in good agreement with theory. The formation of toroids 



on the surface appears to be a first order transition as evidenced by the bimodal 
distribution in n(x). The theoretical result underestimates the simulated n{x) for 
x <C b and follows n(x) closely for x > b; the density calculated exactly agrees well 
with n(x) in the range x <C b. Chain-surface interaction is probed by subjecting 
the surface structures to a pulling force /. The average extension (x(f)) as a func- 
tion of / exhibit sigmoidal profile with sharp all-or-none transition at the unfolding 
force threshold f=f c which increases for more structured states. Simulated 
compare well with the theoretical predictions. The critical force f c is a function of 
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l s /l c for a fixed temperature, where l c and l s are the length scales that express the 
strength of the intramolecular and SC-surface attraction, respectively. For a fixed 
l s , f c increases as l p decreases. 



I. INTRODUCTION 



Interactions between biomolecules and surfaces are important in a number of biological phe- 
nomena. Binding and unbinding of proteins from macromolecular complexes are involved in 
the regulation of biological functions j]} 0, [J. Adsorption of fibrinogen influences the adhe- 
sion of leukocytes, microphages or platelets. In addition, interaction between proteins, DNA 
and RNA are mediated by biological membranes jj, 0, 0| . In the crowded cellular conditions, 
protein-protein and DNA-protein interactions take place in confined geometries in which surface 
interactions are vital. For instance, interaction between P-selectin receptors and their specific 
ligands is mediated by a flat and shallow binding interface [?], Isj]. Besides these situations, 
which are obviously relevant in biology, there are a number of situations in polymer science 
where interactions with surfaces are important 0, Q, 11,^^. These include nanolubrication 
that involve interaction between surfaces that are mediated by polymers. Design of nanoscale 
materials and biologically inspired self-assembling systems also requires an understanding of 
how heteropolymers and biomolecules interact with surfaces. Recent advances in atomic force 
microscopy [l, 13, Q| has allowed a direct probe of the energetics of interaction between ad- 
sorbed proteins with other biomolecules , 19]. The potential applications of 
polymer-surface models to a number of problems has prompted us to develop a theoretical ap- 
proach which can be used in conjunction with AFM experiments to decipher biomolecule-surface 
interactions. 

There have been numerous studies of adsorption of flexible polymers adsorbed on solid sur- 



20]. 



faces which find applications in many aspects of colloidal and interface science 11, I2I 
However, many biomolecules, includin g D NA, RNA, and proteins, are better described using 
worm-like chain (WLC) models Q, 21. Thus, it is important to provide a theoretical de- 



scription of the interactions between semiflexible chains [22J, |23j, |2J] and interfaces. The purpose 
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of this paper is to address the following specific questions: (i) It is known that DNA, a semi- 
flexible poly mer, undergoes a coil-globule transition in the presence of osmolytes or multivalent 



cations 



25| . Simulations of semiflexible chains in poor solvents [20J, |26| have been used to 



understand the kinetics and pathways of transitions from extended conformations to collapsed 
toroidal structures. The coil-globule transition in stiff chains in the bulk occurs through a series 



of metastable racquet structures 
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261 ] . How does the interaction with the surface alter the 



morphology and kinetics of such transitions? This question is relevant even for DNA collapse 
in cells where the DNA compaction takes place in the presence of interactions with their large 
biomolecules in restricted spaces; (ii) AFM experiments are likely to provide the most direct 
data for the strength of interaction between semiflexible biomolecules. In these experiments one 
of the molecules of interest is anchored onto the surface while force is applied to the end of the 
other. The unbinding force can be calculated from the force-extension profiles. These exper- 
iments raise the question, namely, what are the adhesive forces between semiflexible polymer 
and a surface? We address these questions using theory 
model interacting with a solid surface. 
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2811 and simulations for a WLC 



In the absence of the surface the morphologies of semiflexible chain (SC) is determined by 
thermal fluctuations and an interplay of the chain persistence length l p and intramolecular 
condensation length l c =\Jl P kBT /u m where T is the temperature and u m is the effective in- 
tramolecular attractive energy per unit length j^J. In the presence of a surface another length 
scale l s =^/lpkBT/u a ds, where u a( is is the attractive SC-surface interaction energy per unit length, 
plays an essential role in the determination of the structures. The interplay of l p , l c , and l s will 
determine the morphology of the surface-induced structures. It also follows that the response 
to applied force measured in terms of force-extension profiles will depend on l p , l c , and l s . In 
this paper we explore a range of values of fc#T and l s to predict the force-extension curves for 
semiflexible chains in poor solvent. 



II. THEORY 



Consider a semiflexible chain interacting with a flat surface with the SC-surface potential 
being U a d s - Force i=f-n x is applied to one end of the chain (Fig. 1). Equilibrium chain 
configuration is described using the conditional probability G(xjv, f) of finding the tagged N th 
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monomer at x^r given that the monomer xi is anchored at the surface, where x=(r, n) includes 
position vector r=(x,y,z) and orientation vector n, respectively. Due to axial symmetry, the 
free end orientation is specified by the angle between its tangent vector and the x-axis and 
distance from the surface x (Fig. 1), and the conditional probability G(xn, On] #ii /) can be 
used instead of G(xjv, xi; f). In the limit L-^oo, G(xn, On', X\, 9\\ f) is dominated by the ground 
state ipo, so that 

G(x N , N ; xi, Or, f) « VoOi, Oi, f)il>l(x N , N ; f) exp [-/3JVe ] (1) 

€q=Eq/N is the equilibrium free energy per monomer and (5=1/ ksT. If the SC is modeled as 
a worm-like chain (WLC), then if) Q satisfies (22 
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- 2/ p7 ^ + (1 - 7 2 )^ - 2 7 ^ + /5/ p7 /Vo = P(u ads (x) - e )^o (2) 

where u a d s = U a d s /N is the adsorption potential per monomer and 7 = cos [0] (see Fig. 1). To 
mimic the Lennard- Jones chain-surface attractive interaction used in the Langevin simulations 
(see Section III), we employ a piece- wise continuous potential, i.e. u a d s =oo for x<0, u a d s =— A 
for 0<x<b and u a d s =0 for x>b. The monomer density of the adsorbed structures in the absence 
of force, 

n(x) = [ ddifi${x,6) (3) 



normalized as J dxn(x)=N, is calculated by solving Eq. (J2J) without the last term l3l p ^fipQ. 
The perturbative solution of Eq. (J2J) in the absence of the flippy f term, due to Kuznetsov and 



Sung 



27j . to the first order in correlation length parameter T)=(4:lp/ipo)\cP ip /dx | is outlined in 



Appendix A. The solution is 



M*) = Kb - x){C x sin (^P-) + C 2 sin ) (4) 



^o(^) = Coh(x — 6) exp 



x 



2/ p 



where /i(x) is the Heaviside function, and Co, Ci, C2 and C3 are constant coefficients; m 1; m 2 
and k are given by 

mi, 2 = ^(2 + in )(6 + m ) I 1 T + I {.',} 
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k = y(2 + </w)(6 + oui ) 1-1 + 

and 4>in=/3(u a ds(x)— eo)<0, 4>out=— e o>0- By using two continuity requirements (jA9|) and the 
normalization, we can obtain, respectively, eo and Co and one of Ci, or C3. However, the 
two free constants are to be chosen such that a minimum of eo is obtained. The minimal free 
energy corresponding to the ground state for x < b is attained for (i) Ci^O, C 2 =C 3 =0 (i.e. the 
state with m—mi) and (ii) 6*27^0, Ci=C3=0 (171=1712). 

Perturbative solution of Eq. (J2J) ignores variation of ipo on 9. Indeed, when x3>/ p , ip becomes 
nearly isotropic, ip(x, 9)=if)(x). However, when x~o<^Z p , ip should strongly depend on the 
angle G = it/2 + 8 between free end of the chain and the surface (Fig. 1). In this range 
—jd/dx^-Qd/dx, (1 — 7 2 ) d 2 /d'y 2 ^d 2 /dQ 2 , 7<9 /cfy^O, 7/^—6/ and Eq. (J2J) simplifies, i.e. 

dipo . 1 d 2 ij 1 f , P , , n x , / r n 



The methodology for solving Eq. (JOJ) has been presented by Semenov in Ref. j28] and is outlined 
in Appendix B. The general solution for x>b (u ads =0) is 

Mx,e)=^^Ll,^) ( 7) 

71=0,1 ^ ' 

where Co, C\ are constants and the confluent hypergeometric function ty(X)W,z) is 
p(y, q>, z)=Y^ dr r x_1 (l + r) w_x ~ 1 e" TZ where T(x)= J °°(irr x ~ 1 exp [— r] is Gamma function 
[29| . To describe the chain in the range x < b, we assume that ipo is of the form (J7J) and Co and 
Ci depend on x. Substituting Eq. (jZJ) into Eq. (JHJ), we obtain: 



/2/ \ 1/3 

0(z) MJ (F 00 C + F^ 1 ^) (8) 



(iCo , / s f 2/^ 
dx 



dC\ . . / 2/p 



1/3 

— ) (Fio^Co + Fnd) 
ax \ x J 

where F nm =(g m /n, f n )/{f n , /„), n,m = 0,1, and (#„, # m ) = j^dnne'^ 19 g n (K) g m (K) (see Eq. 
(|B2|) ). We solve Eqs. (JHJ) subject to the condition 0(x)=0j n for x<6 and 0(x)=</w for x>6. 
From the solutions of Eqs. (JHJ) in Appendix B we obtain: 

C{ n (x) = c 2 e^ J (-3y/=Dx^ 1 $ (p + |, |, -3>/=£a;§) - ar) (9) 
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where y/—D = 4> ou t{2l p ) 3 ^F m Fi$ — F n F 00 and h(x) is Heaviside step function. In Eq. (JHJ) the 
Kummer function is defined by /, x) = l+J2m=i^) 11 ' ^7 w hh (k)o=l, {k)\=k and 

{k) m =k(k + l)...(k + m-l) . We get C (x) by substituting Eqs. © for C{ n and Cf* into 
the second Eq. (jHJ) and ^o(^) ©j /) can now be obtained by using Eq. (JJJ)). 

In the presence of pulling force, -00 is nearly isotropic, i.e. ipo{ x , Q] /)~ V'oC^i /)• This allows us 
to analyze force-extension profiles by employing perturbative treatment outlined above. Solution 
to Eq. (j2J) is given by 

M*;f) = M*)e* 0fx (10) 

where ipo(x) is given by Eqs. (@J). The average extension as a function of applied force can be 
computed using 

w» - bk>l z(/) < n > 

where the partition function Z(f) is Z(/)= f d6>jv / J d x Nj dxiG(xN, 9 N ; x±, 9\\ f). 

The perturbation theory is strictly valid only when the condensation length / c 3>/ p . In practice 
we find that the first order perturbation theory gives results that are in very good agreement 
with simulations even when Z c ~/ P . Kuznetsov and Sung also discovered that the perturbation 
theory is remarkably successful outside the regime of applicability j^ . 

III. LANGEVIN DYNAMICS SIMULATIONS 

We model a semiflexible chain (SC) by iV=100 connected beads of bond length a and the 
contour length L=100a. In the absence of U a ds and /=0, we assume that the dynamics is 
governed by the overdamped Langevin equation: 

4*=-S + * (t) (12) 

where £ is the friction coefficient, U=U c h a in=Ubond+Ubend+ULj is chain internal energy due to 
bond potential Ub on d, bend potential Ubend and interbead interaction potential Ulj (hydrody- 
namic interactions are ignored). The random force gj{t) obeys Gaussian statistics, 

( gj (t)) = 0, (gi(t) Sj (t')) = 6k B T^5(t - t') (13) 
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We solve Eq. (|T2j) for each Xj with unit tangent vector Uj=(x.j +1 — Xy)/a, where j=l, 2, . . . ,N. 
The stretching potential t4 on d is 

AT-l 

2<X 2 



^ond = ^^(|x J -X J+1 | 2 - ( T) 2 , (14) 

where A and a are constants, and 

£4end = - ^ (1 + COS [V2i,j+l]) 2 (15) 
3=1 

where the constant S is a measure of chain stiffhes, and cos [ l Pj,j+i\= ( x j+i — x j)( x j-i — x j)/°" 2 is 
the bend angle. The interaction between beads is given by the 12-6 Lennard- Jones potential, 



12 

c \ _ / a 



AXjj / \ AXjj 



(16) 

\ ZAX;„- / \ iAXii / 

where Ax^- is the distance between beads i and j, and 5 is the magnitude of interaction. Ulj is 
an effective interaction that accounts for excluded volume interactions and counterion induced 
attraction which in DNA is due to screening of the charges. The persistence length of the chain l p 
can be roughly estimated by using l p =a/ (1 — cos [(<fj,j+i)]) where (Vj,j+i) = (N —l)^ 1 J2j ! =i i Pj,j+i 
is the average angle between adjacent beeds. 

Similar models have been used in previous studies to probe the chain collapse in poor solvents 
j3, In the presence of the adsorbing surface the motion of j-th bead is governed by Eq. 
(fT2"|) with U=U c hain+U a ds, where U a ds is the surface-SC potential, 

In Eq. (JT7J) Axj is the bead-surface distance and A and b are, respectively, the depth and range of 
the attractive forces. We set B = 1.0, a=a=l and b=3a, and use A=400B, S=30B, 60S, 120.B 
and A=1.5.B, 2.05, 2.55. This makes Ulj, Ubond, Ub en d and U a ds to scale in units of eh=ksT 
and ei=a is the unit length. The choice A=400B allows for 5 percent of thermal fluctuations in 
the bond distance and permits us to run simulations with longer time steps without affecting 
bond relaxation time. The unit of time is r=^cx 2 /e/ l , where £= 44.0 is the friction coefficient of 
the chain in water at T = 3007^. The system of equations (J 12)) is integrated with a step size 
5t=2xl0~ 2 r and the total time is t=N to tSt where N to t is the number of integration steps. We 
express time either in units of r or in terms of N tot . 
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IV. RESULTS 



Surface-induced structural transitions: 



It is known that in the absence of the surface SC undergoes a collapse transition when the 
solvent is poor i.e., when the attractive monomer-monomer interactions dominate (Eq. (|16p) 
so that l c >l p . The collapse is a result of a competition between intramolecular attraction and 
bending energy due to chain stiffness. Unlike in flexible polymers, the low energy collapsed 
conformation is a torus which maximizes intramolecular contacts and minimizes the bending 
penalty. Before simulating the force-extension curves of adsorbed SC it is necessary to charac- 
terize the structures that are obtained when interacting with the surface. 

To simulate the low free energy structures that result in the presence of the surface, we 
first thermalized an extended chain at high temperature fc#T=3.0 for iV tot =lxl0 6 steps. By 
gradually decreasing the temperature bulk structures were thermalized for (1 — 10)xl0 7 steps 
and used in adsorption experiments. Interactions with the attractive surface was switched on at 
distance Ax=2b away from the bead with shortest x and the SC was adsorbed onto the surface 
one bead at a time. The structures were allowed to relax for ~1— 20xl0 6 steps depending on 
ksT, S and A. Progress of adsorption was monitored by analyzing time traces of Ulj, U a d s , 
U, and the radius of gyration R g of the SC. We generated 500 adsorbed structures at k B T=1.0, 
1.25 and 1.5 for S=30B, 605, 1205 and A=1.55, 2.05 and 2.55. 

Typical structures are presented in Fig. 2. Geometry of the SC adsorbed onto the surface 
ranges from partially or fully extended configuration with l p /a ~ 18 to partially structured 
one-, two- and three-racquet states with Z p /a~16.0, 15.5 and 15.0, respectively, to fully ordered 
toroidal states with l p /a ~ 13.5. Similar structures have been observed in recent studies of 



collapse of semiflexible chains in the bulk 
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261 ] . For the interaction parameters used in our 



simulations i s /i c ~o(l). Thus, the attractive SC-surface interaction facilitates adsorption of the 
SC without significantly altering its morphology compared to the bulk case. For l s ~^>l c the lowest 
free energy structures are extended. 

To compare the kinetics of structure formation on the surface and in the bulk we also sim- 
ulated collapsed structures in the absence of the adsorbing surface. By analyzing the temporal 
profiles of R g , Ulj and l p , we found that on average, chains attain structured configurations on 
a faster timescale when adsorbed on the surface. The search for the ground (toroidal) state is 
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more efficient when the chain is constrained to evolve on the two-dimensional surface where the 
SC quickly minimizes its free energy in reduced d=2-space by sliding surface motion (lateral 
diffusion) . 

B. Kinetics of surface-induced ordering: 

Typically, surface-induced ordered structures form by a two step process B -^So—>-St- Start- 
ing from the bulk state B , extended surface transient S emerges during the fast first step 
with the B ^So transition occuring within N to t=l — 3xl0 6 . In the slower second step So^S t , 
extended transient structures explore the free energy landscape in search of the toroidal state St 
which occurs in about 3 — 20xl0 6 steps depending on S, A and temperature. Transition from 
So to S t is realized via rapid formation of either a surface loop or an intermediate toroid-like 
motif with larger R g (smaller winding number) or through a sequence of longer lived racquet 
states Sq^Si^. . .— ^St, where S n , n=0,l,2,. . . denotes conformations with number of racquets 
equal to zero (extended chain) one, two, etc. 

The number of "metastable" racquets depends on chain flexibility. We observed configu- 
rations with n=6 for S=30B and ksT '=1.25. Simulated profiles of R g , Ulj, U a d s and U in- 
dicate that evolution from extended to toroidal states follows several pathways. Four out of 
five simulation runs followed the scheme outlined above. Similar diverse pathways have been 
observed by Noguchi and Yoshikawa who recorded the lifetime of intermediates species for 
about iV to4 =2.0xi0 5 . Our results indicate that attractive surface forces increase the lifetimes 
of metastable intermediates for stiff chains at low temperature. In few simulation runs toroidal 
structures were not observed during as many as 20xl0 6 steps. Hence, attractive surface forces 
facilitate formation of toroidal state primarily when formation of toroid-like intermediate motif 
is involved. 

The dynamics of R g , Ulj, U a d s and U for the structures in Fig. 2 show that increasingly 
more ordered states are also energetically favorable (Fig. 3). R g , U and Ulj decrease and U a bs 
increases in the sequence Sq— hSi— hS-2— >S^— >St- R g fluctuates around larger values for extended 
states. Variations in Ulj, U a d s and U increasing in the sequence of Si—tSz—tS^—tSt transitions 
are due to formation of SC-surface contacts. For the structures in Fig. 2, the formation of Si 
at N tot ^2.0xl0 6 is mediated by a surface-loop motif followed by slow sliding motion; 5*2 forms 
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early at iV to ^5.0xl0 5 and remains unchanged (time dependence of R g or Ulj). The dynamics 
of Ulj, U a ds and U show formation of S3 via 5*1 at A^ ot ?»1.0xl0 5 followed by transition S\-^Sz 
at N tot ~7.0 x 10 6 . Similarly, traces of same quantities for St point at three step transition, 
So—>Si—>Sa—>S t occuring respectively at iV toi ^5.0x 10 5 , 1.5xl0 6 and 4.0xl0 6 , followed by chain 
compaction due sliding motion. 



In agreement with theoretical arguments 



2 a ] monomer profiles of stiff chains (S=120B, 



l p /b^$> 1) are described by the succession of short near-surface loops of length <C l p /a between 
chain-surface contact and by the combination of short and long loops of the length ^> l p / a for 
S=30B and Z p /fe~l. Decrease in l p /b and temperature favors the formation of chain-surface 
contacts by enabling more beads to be inside the range of surface forces. This results in the 
formation of higher ordered states S4, S$, Sq and S t . In contrast, at higher temperatures and 
increased A and S, surface structures with increased conformational free energy become unstable 
and unfold into extended configurations (data are not shown). We quantified the geometry of 
surface structures (Fig. 2) by binning bead-surface distances Xj into the density histogram 
n(x). The monomer density profiles for A=1.5B and kf,T =1.0 are compared in Fig. 4 for 
S=70B (left) and S=50B (right panels). Transition from less structured to more structured 
states is accompanied by an increased ratio of the number of bead-bead to bead-surface contacts. 
The density distribution n(x) is single-peaked at x=b/2 and decays to zero as x— >b for extended 
states and increases its density at a;~6 in the sequence Si—*S2—*S3—>St. 



C. Forced unfolding of surface adsorbed structures: 

To unfold the surface-ordered structures, these structures were initially allowed to thermalize 
at /cbT=1.0 for N to t=2 x 10 6 . We then ancored the C-terminus of the chain at the surface and 
pulled its iV-terminus with constant force / via the harmonic spring with the spring constant 
fc sp =0.36pA^nm~ 1 in the direction perpendicular to the surface. Simulation runs were terminated 
after evolution of chain extension x(N to t) had reached equilibrium. x(N to t) of the structures of 
Fig. 2 are presented in Fig. 5 for f=9.75pN , 18.3piV, 24ApN and 30.5pN. Chain extension 
reaches saturation plateau in the first 8xl0 7 steps as the chain restoring force approaches /. 
Not unexpectedly, the unfolding threshold force increases as the extent of ordering decreases in 
the sequence So^Si^S 2 ^S3^S t . At f=9.75pN only S unfolds in l.OxlO 7 steps. When the 
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force is increased to f=18.3pN, Sq, Si, 5 2 and 5 3 unbind from the surface in 3.5xl0 7 , 3.6xl0 7 , 
4.0xl0 7 and 6.0xl0 7 steps, respectively. At f=2AApN all structures reach the stretched state 
in 2— 4xl0 7 steps. From the dynamical trajectories of x obtained for A=1.55 and fceT=1.0, we 
constructed the average extension (x) as a function of /. In Fig. 6 we compare (x) vs / traces 
for extended, one-racquet, three-racquet and toroidal structures of Fig. 2 (5=1205, top panel) 
and more flexible four-, five- and seven-racquet and toroidal conformations obtained for 5=305. 
Unbinding of surface- anchored structures undergo a highly cooperative all-or-none transition as 
the unfolding force threshold f=f c is increased from 7.3— 15.8pN (5=1205) to 15.9— 17.7pN 
(5=305) for more compact racquet and toroidal states. 

D. Comparison between theory and simulations: 

We analyzed the simulation results for the monomer density and the averaged extension 
as a function of the pulling force by using perturbative treatment (see Eqs. ©-©) i n the 
entire range of x<L. For the proximal limit x^l p <L we use the exact expressions in Eqs. (JTj) 
and (jHJ). Density distributions and force-extension profiles for the extended conformation were 
approximated by chosing the ground state with m=m\ (mi<m2, see Eq. ©). The choice 
m=mi corresponds to isotropic-like unstructured surface state with no prefered orientation of 
the chain beads. Histograms of structured two-, three-racquet and toroidal conformations were 
analyzed with the choice m=m 2 corresponding to nematic-like ordered states |27(. To account 
for the difference between the shape of attractive potential U a d s used in the simulations and the 
theoretical calculation we used, in the actual fit, the rescaled potential depth Ar=rA S j m for the 
same range bT=b S i m =l, where r=(5A*r) _1 f^°dxu a d s (x) is the ratio between volume of Lennard- 
Jones attractive layer and bAj- used in theory. The density profiles n(x) for known values of 
b, ksT and Ay were fitted to the simulated monomer density histograms and force-extension 
profiles to obtain parameters eo (Eq. (£Q)) and l p . The theoretical results for the density n{x) 
and the average extension (x(f)) computed from Eqs. Q and (JUJ), respectively, using these 
parameters are shown in Figs. 4 and 6. 

Monomer density distributions: Although the theoretical results for n(x) slightly underes- 
timate the simulated density for structured states for b<x/a<1.5 and underestimates it for 
x/a>1.6, the agreement between perturbation theory and simulation data is surprisingly good 
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in the range of x/a<b (Fig. 4). The agreement between theory and simulations improves 
for more structured racquet and toroidal conformations. In particular, the theoretical profiles 
capture the positions of density peaks both inside the layer at x/a^0.5 and at the boundary 
x/a~b. Although there is some residual density at large x/a due to thermal fluctuations of chain 
ends, especially for less structured extended and racquet configurations, the ground state domi- 
nance approximation is clearly valid. The theoretically estimated conformational free energy per 
monomer and persistence length for structures S ^S 2 ^S 3 ^S t of Fig. 4 decrease respectively 
as eo/fcBrw-0.21->-0.23->-0.24->-0.25 and Zp/awll.7-^11.2-^10.4-^10.2 (for k B T=l.O, 
left panels), and eo/JfeBr«-0.23->-0.25->-0.26->-0.27 and Z„/a«11.3->11.0->10.2->10.1 (for 
k B T =1.25, right panels). Not surprisingly, both eo and l p decrease for the same structures 
as k B T is increased because of enhanced chain flexibility. In the proximal region, the exact 
calculation of n(s) for 0<x/a<0.5 for the same structure sequence shows a better agreement 
with the simulated results. The fit parameters are e /k B Tm— 0.2^— 0.24^— 0.25^— 0.28 and 
Zp/a«12.3->11.8->11.0->10.8 for k B T=1.0, and e /fc B T«-0.19-»-0.22-»-0.24-»-0.26 and 
Zp/a«12.0->11.6->10.8->10.6 for fc B T=1.25. 

Force- extension curves: Apart from small deviations around the unfolding threshold 
forces for all simulated surface structures, the fit of theoretical curves of the average ex- 
tension vs pulling force to simulated data points shows excellent agreement between the- 
ory and simulations. The theoretical (x(f)) curves calculated using perturbation theory 
follow closely the simulated force-extension profiles both for 5=1 205 and S=30B espe- 
cially below (x/L<0.1) and above (x/L>0.9). The unbinding threshold forces increase as 
7.5pN<10.5pN<12.5pN<16.5pN in the sequence S ^Si^S 3 ^S t (S=120B, top panel in Fig. 
6) and as 15piV<15.5piV<16.5piV<17.5piV in the sequence S4— >55— >5V— >£j (S—30B, bottom 
panel in Fig. 6). This implies that more flexible and/or more structured surface chains are 
harder to unfold. However, "all-or-none" type of simulated unfolding transition shows sharper 
growth than predicted by the theory. The theoretically estimated conformational free energy 
per monomer and persistence length for structures So^Si^S^St decrease respectively as 
eoAijT^-0. 12^-0. 17^-0.20^-0.22 and Z p /a^l5^14.5^14.25^12.1 (top panel). For the 
structures S4— ><S5— ><SV— e decreases as e /k B T^i— 0.134^— 0.136^— 0.141— >— 0.148 and 
Z p /a^8.2-^8.1^8.0^7.9 (bottom panel). Here too, increased chain flexibility decreases l p and 
lowers e . 
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V. CONCLUSIONS 



To provide insights into interactions between biomolecules interacting with membranes we 
have considered collapse and forced-unbinding of semiflexible chains (SC) in the presence of an 
adsorbing surface. The interaction of SC modeled using WLC, which describes well many of the 
physical properties of DNA Q, RNA Q, and proteins [3JJ], with a surface into which the SC 
can adsorb, is studied using theory and simulations. The morphologies of the SC in the presence 
of an adsorbing potential is described in terms of three length scales, namely, l p , l s , and l c . By 
restricting ourselves to i c ~Z s we have studied the effect of interaction with the surface on coil- 
toroidal transition in DNA like chains. The simulations show that the rate of toroid formation is 
impeded compared to the bulk because interaction with the surface stabilizes many metastable 
racquet-like structures (Fig. 1). The simulated equilibrium density profiles show that as the 
range of surface-SC interaction increases and temperature decreases, which leads to a decrease 
in lp/l s , ordered structures form. The peak of n{x) at x~b (the range of interaction) grows as l p 
decreases. The bimodality in the n{x) distribution function suggests that the surface-induced 
toroid formation is a first order transition. The perturbative calculation reproduces qualitatively 
all the features in the simulated density profiles. 

We also considered the peeling and unbinding of adsorbed structures by applying force. 
These results, which are of direct relevance to AFM experiments [32], show that the forced- 
unbinding transition is surprisingly highly cooperative. For all structures (racquet-like and 
toroids) unbinding occurs over a narrow force range. The magnitude of the critical force f c 
for a fixed value of T and l s increases as l p decreases. From general considerations we expect 
that f c should be described by a scaling a function g(y) where y=l s /l p for a fixed T. When 
y<y c (a critical value), then adsorbtion is not free energetically favored. When y>y c , then 
f c should increase by an increasing function of y. The increase in f c can be achieved either 
by increasing l s for a fixed l p or by decreasing l p for fixed l s . Additional work is required to 
elucidate the nature of the scaling function g(y). Quite surprisingly, we find that the force- 
extension profiles can be calculated by using a simple perturbation theory even though the 
nature of the unbinding transition is abrupt. The present work shows that global properties of 
force-extension characteristics of adsorbed biomolecules can be nearly quantitatively predicted 
using the proposed theory. 
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It is now well established that elastic response of DNA, in the absence of interaction with 
surfaces, depends sensitively on the nature and concentration of counterions 



32 



331 ] . Our work 



shows that the force-extension curves in the presence of a surface to which DNA is bound depends 
not only on l s but also on the morphology of the adsorbed structures. The novel prediction that 
forced unbinding should occur cooperatively by a first-order phase transition can be probed 
using single molecule experiments. 
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APPENDIX A: PERTURB ATIVE TREATMENT OF ADSORBED CHAIN 

STATISTICS 

We expand ipo (Eq. ([T])) in terms of the Legendre polynomials, i.e. 

oo 

VoM) = X>^) p *( cos [0D ( A1 ) 

By using the following equations, 

±P M = -i {i + mh) and P 1 ( 1 )P i ( 1 )= {i + 1)P 'it! iPi - lM < A2 > 

Gt7 ll + 1 

we transform Eq. (J2J) without term pi p jfip into {i > 0) 

+ (3(u ads - eg) . . i dj)j-i(x) _ i + 1 dip i+1 (x) 

2l p [X) 2i - 1 dx 2i + 3 dx [ ' 

To the first few orders we have: 

, = 21 p d ^ (a = o) (A4) 
° ?>l3{uads - eo) dx 

^ 2/ p dipo 4/p <hp2 = 

2 + (3(u ads - e ) dx 5(2 + (3(u ads - e )) dx 

^ 4/p (hpi 6lp di/js U-2) 

3(6 + (3(u ads - e )) dx 7(6 + (3{u ads - e )) dx 
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ipi is given by the second Eq. (jA4|) with the second term determined by -02 which is of the order 
of 8lpd 3 tp /dx 3 . Neglecting this order for 7?<Cl we obtain: 



2L 



dipo 



(A5) 



2 + (3(u ads - e ) dx 

Including the second term in the third Eq. (jA4|) and using Eq. (jA5|) we obtain the first-order 
perturbation equation for ipi, 



2l p dipo 



AL 



d 



Alp d 



21 p dipo 
3(6 + 0) dx \2 + <ix 



(A6) 



2 + (f) dx 5(2 + 0) ofe 
where = f3{u a d s — eo)- Substitute Eq. (jA6|) into the first Eq. (jA4|) we arrive at the first-order 
perturbative equation for ip : 

d 2 i/j 



64/J d^b , All d\« 



(A7) 



450(2 + 0) 2 (6 + 0) dx 4 30(2 + 0) dx 2 
For the class of potentials considered here, the physical solution of Eq. (jA7|) that satisfies the 
boundary conditions, 

d" 



ip (x = 0) = 0, ip (x -> oo) 
continuity requirements at x — b, 



0. 



dx' 



1 x— »oo 



d 



0, n=l,2, 



d 



ipo(x -»• & - 0) = ^o(a; -> & + 0), — ^oU-^6-o = -j-^oU-^b+o 



dx 



<ix 



(A9) 



and normalization condition is given by Eq. (0J). 



APPENDIX B: EXACT TREATMENT OF THE CHAIN DISTRIBUTION IN THE 

PROXIMAL RANGE 

Let us first consider the non-adsorbed chain in the presence of weak potential u a d s — >0- As- 
suming a self-similar distribution, i/)q(x,Q) =x a g(n) where k=0(2/ p /x) 1//3 , we rewrite Eq. 
with M a ds=eo=0 as an eigenvalue problem for #(k), 

-\^l + 1^L = a (Bl) 
k 8k 2 3 Ok 

Upon substitution z=k 3 /9, Eq. (jBl|) reduces to the following equation: 

d 2 g (2 \ da , s 

3- 2 h! +as = (B2) 
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Under condition g(z—>-oo)^0 the only solution to Eq. (jB2|) is g(z)~^(— a, 2/3, z) where 
^{Xi u i z ) an d T(x) are defined in the main text. ip(x— >0,6<0) — >0 defines the spectrum of 
eigenvalues a n =l/6— n, where n = 0,±1,±2,. . . (see Appendix B in Ref. j^)- The requirement 
that ip does not have knots is satisfied for n = (a=l/6) and n—1 (a—— 5/6), and the general 
solution for u a d s =0 in the region x>b is given by Eq. (jZJ) of the main text. 

To solve Eqs. (jHJ we substitute C from the second equation into the first equation and 
multiply by x 2//3 . We obtain: 

xA/z <m _ 2 usdCi + + Dx2/ = o 

ax z 3 ax 

where .D = Foo-fii — -F 01 F 10 and F nm = F nm (f)(2l p ) 1 ' 3 , n, m = 0, 1. Substituting y = x 1//3 into 
Eq. (|B3j) and multiplying it by y 2 , we get: 

_ dc± + - u + Ci = Q 

Using z = y 2 allows us to rewrite Eq. (jB4|) as 

l2Z ^i + + (7oZ + A)cfl = (B5) 



where 7o=9.D, p =9Fu, (3\=—A and 72=4. The general solution of Eq. (|B5|I is given by 



= e f^* s [ Cl $(p, -±, - 3v ^Dxi) + c 2 (-3V=Pxt)f $(p + ^, ^, -Sv^xi)] (B6) 

where Ci, c 2 are constants and p=(3Fu—2\/—D)/4\/—D. In Eq. f)B6|) x) is the Kummer 

series defined in the text. In the range 0<x<o, ipi(x,Q=0) diverges as x^O. To avoid this 
divergence we require that Ci(x=0)=0. This is satisfied when ci=0. To insure that Ci(x)— >0 as 
x^oo for x>6, we set c 2 =0. Substituting <p—(pi n and (p=4> ou t into solutions for 0<x<6 and x>b 
and using formulas I, x)=e x &(l — k,l,x), ^ r $(/c, Z, x)=-^rp-$(/c + m, / +m,x) we obtain 
Eqs. ©• 



_> 



1 „ / 2 . . „ , — 2 . 3 T . 3 5 



[1] Liphardt, G.; Smith, D.; Bustamante, C. Curr. Opin. Struct. Biol. 2000, 10, 279; Allemand, J.-F.; 

Bensimon, D.; Croquette, V. ibid 2003, 13, 266 
[2] Yin, H.; Wang, M. D.; Svoboda, K.; Landick, R.; Block, S. M.; Gelles, J. Science 1995, 270, 1653; 

Bianco, P. R.; Brewer, L. R.; Corzett, M.; Balhorn, R.; Yeh, Y.; Kowalczykowski, S. C; Baskin, 

R. J. Nature 2001, 409, 374 



17 



[3] Alberts, B. Molecular Biology of the Cell, 1st ed.; Garland Publishing: New York, 1994 

[4] Labeit, S.; Kolmerer, B. Science 1995, 270, 293; Minajeva, A.; Kulke, M.; Fernandez, J. M.; Linke, 

W. A. Biophys. J. 2001, 80, 1442 
[5] Henrickson, S. E.; Misakian, M.; Robertson B.; Kasianowicz, J. J. Phys. Rev. Lett. 2000, 85, 3057; 

Kasianowicz, J. J.; Brandin, E.; Branton, D.; Deamer, D. W. Proc. Natl. Acad. Sci. USA 1996, 

93, 13770 

[6] Sung, W.; Park, P. J. Phys. Rev. Lett. 1996, 77, 783; Muthukumar, M. J. Chem. Phys. 2003, 118, 
5174; 1999, 111, 10371 

[7] Marshall, B. T.; Long, M.; Piper, J. W.; Yago, T.; McEever, R. P.; Zhu, C. Nature 2003, 423, 
190; Ramachandran, V.; Yago, T.; Epperson, T. K.; Kobzdej, M. M.; Nollert, M. U.; Cummings, 
R. D.; Zhu, C; McEver, R. P. Proc. Natl. Acad. Sci. USA 2001, 98, 10166 

[8] Leppaenen, A.; White, S. P.; Helin, J.; McEver R. P.; Cummings, R. D. J. Biol. Chem. 2000, 275, 
39569; Somers, W. S.; Tang, J.; Shaw G. D.; Camphausen, R. T. Cell 2000, 103, 467 

[9] de Gennes, P.- J. Scaling Concepts in Polymer Physics, 1st ed.; Cornell University Press: Ithaca, 
1979; Freer, G.; Cohen-Stuart, M.; Scheutjens, J.; Cosgrove, T.; Vincent, B. Polymers at Inter- 
faces, 1st ed.; Chapman and Hall: London, 1993 
[10] Muthukumar, M. J. Chem. Phys. 1987, 86, 7230; 2004, 120, 9343; Wiegel, F. W. J. Phys. A. Math. 
Gen. 1977, 10, 299 

[11] Ellis, M.; Kong C. Y.; Muthukumar, M. J. Chem. Phys. 2000, 112, 8723; McNamara, J.; Kong C. 

Y.; Muthukumar, M. ibid 2002, 117, 5354 
[12] Semenov A. N.; Joanny, J. F. Europhys. Lett. 1995, 29, 279; Semenov, A. N.; Bonet-Avalos, J.; 

Johner, A.; Joanny, J. F. Macromolecules 1996, 29, 2179 
[13] Weisel, J. W.; Shuman, H.; Litvinov, R. I. Curr. Opin. Struct. Biol. 2003, 13, 227; Liphardt, J.; 

Dumont, S.; Smith, S. B.; Tinoco, I.; Bustamante, C. Science 2002, 296, 1832 
[14] Chang, K. C; Tees D. F.; Hammer, D. A. Proc. Natl. Acad. Sci. USA 2000, 97, 11262; Bustanji, 

Y.; Arciola, R. C; Conti, M.; Mandello, E.; Montanaro, L.; Samori, B. ibid 2003, 100, 13292 
[15] Rief, M.; Clausen-Shaunman, H.; Gaub, H. E. Nat. Struct. Biol. 1999, 6, 346; Rief, M.; Gautel, 

M.; Schemmel, A.; Gaub, H. E. Biophys. J. 1998, 75, 3008; 
[16] Ha, T.; Zhuang, X.; Kim, H. D.; Orr, J. W.; Williamson J. R.; Chu, S. Proc. Natl. Acad. Sci. 

USA 1999, 96, 9077; Zhuang, X.; Bartley, L. E.; Babcock, H. P.; Russell, R.; Ha, T.; Herschlag, 



18 



D.; Chu, S. Science 2000, 288, 2048 
[17] Hyeon, C; Thirumalai, D. Proc. Natl. Acad. Sci. USA 2003, 100, 10249 

[18] Barsegov, V.; Thirumalai, D. Proc. Natl. Acad. Sci. USA 2005, 102, 1835; Barsegov, V.; Thiru- 
malai, D. Phys. Rev. Lett, (in press) 

[19] Bustamante, C; Marko, J. F.; Siggia, E. D.; Smith, S. Science 1994, 265, 1599; Smith, D. E.; Chu, 
S. ibid 1998, 281, 1335; Perkins, Th. T.; Smith, D. E.; Chu, S. ibid 1997, 276, 2016 

[20] Schurr, F.; MacKintosh, F. C; Williams, D. R. M. Europhys. Lett. 2000, 51, 279; Montesi, A.; 
Pasquali M.; MacKintosh, F. C. Phys. Rev. E 2004, 69, 021916 

[21] Bouchiat, C; Mezard, M. Phys. Rev. Lett. 1998, 80, 1556; Marko, J.; Siggia, E. D. Macromolecules 
1995, 28, 8759; Wilhelm, J.; Frey, E. Phys. Rev. Lett. 1996, 77, 2581 

[22] Samuel, J.; Sinha, S. Phys. Rev. E 2002, 66, 050801; Krishna, M. M. G.; Samuel, J.; Sinha, S. J. 
Phys. A 2000, 33, 5965 

[23] Lifshitz, I. M. Zh. Eksp. Teor. Fiz. 1968, 55, 2408; Grosberg, A. Yu. Biofizika 1979, 24, 32 (Bio- 
physics 1979, 24, 30) 
[24] Kierfeld, J; Lipowsky, R. J. Phys. A: Math. Gen. 2005, 38, L155 

[25] Post, C. B.; Zimm, B. H. Biopolymers 1979, 18, 1487; Bloomfield, V. A. ibid 1997, 44, 269 
[26] Noguchi, H; Yoshikawa, K. J. Chem. Phys. 2000, 113, 854; Sakaue, T.; Yoshikawa, K. ibid 2002, 
117, 6323 

[27] Kuznetsov, D. B.; Sung, W. J. Chem. Phys. 1997, 107, 4729; Macromolecules 1998, 31, 2679 
[28] Semenov, A. N. Eur. Phys. J. E 2002, 9, 353 

[29] Abramowitz, M.; Stegun, I. Handbook of Mathematical Functions, 9th ed.; Dover: New York, 
1972. 

[30] Caliskan, G. submitted to Phys. Rev. Lett. 

[31] Schuler, B.; Lipman, E. A.; Steinbach, P. J.; Kumke, M.; Eaton, W. A. Proc. Natl. Acad. Sci. 

USA 2005, 102, 2754 
[32] Murayama, Y.; Sakamaki, Y.; Sano, M. Phys. Rev. Lett. 2003, 90, 018102 
[33] Lee, N.; Thirumalai, D. Eur. Phys. J. B 1999, 12, 599 



19 



FIGURE CAPTIONS 

Figure 1: Schematic of a semiflexible chain (blue) adsorbed on the surface (yellow). The 
free end (r, n) makes an angle 6 l =arccos [n • n x /|n| • In^l] with the direction n x of the pulling 
force f=fn x . For clarity the chain is shown as extended which is realized only when the SC- 
surface interaction is strong. The interaction between the monomers of the chain and the 
surface is attractive in the range 0<a;<6 where x is the distance perpendicular to the surface. 
The strength of the interaction is A. In the Langevin simulations we replace the square well 
potential by the Lennard- Jones potential (Eq. (fTTj) ). 

Figure 2: Top view of the typical structures (blue) adsorbed on the surface (yellow) 
for 5=1205, A=1.5.B and fcsT=1.0. Extended, one-racquet, two-racquet, three-racquet and 
toroidal structures are obtained in a single trajectory that is terminated at t=4xl0 _5 r. The 
equilibrium structure under these conditions is the toroid. 

Figure 3: Dependence of radius of gyration R g /a (top left), intramolecular attractive 
interaction Ulj (Eq. ([Tfi])). surface potential U a d s (Eq. (JTTJ) ) , and the internal energy U 
(=Ubend+Ubond+ULj) displayed as functions of time, measured in units of r. The five curves 
in each panel correspond to extended (black), one-racquet (red), two-racquet (green), three- 
racquet (blue) and toroidal (magenta) structures of Fig. 2. 

Figure 4: The average monomer density profiles n(x) as a function of x/b for A=1.5B and 
/cbT=1.0 for extended, two-racquet, three-racquet and toroidal states. The left panel is for 
S=70B and the results for S=50B are shown on the right. Solid lines and dotted lines represent 
the results obtained using perturbative and exact theory, respectively. 

Figure 5: Dynamics of extension x (in units of a) for a few trajectories at different values 
of / applied to the chain ends of structures shown in Fig. 2. Time t is expressed in units of r. 
The colors correspond to the caption in Fig. 3. The values of / are displayed in the panels. 

Figure 6: The averaged reduced extension (x)/L as a function of constant force / simulated 
for A=1.5.B and ksT— 1.0 for structures in Fig. 2 (S=120B, top). The bottom panel shows 
force-extension profiles for four-, five-, seven-racquet and toroidal configurations obtained for 
5=30.6 (bottom panel). Data points for extended (four-racquet), one-racquet (five-racquet), 
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three-racquet (seven-racquet) and toroidal structures are given by red, green, blue and black 
circles, respectively. Theoretical curves for these structures are given respectively by solid, 
dotted, dashed and dot-dashed lines. 
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Ucds = 0, x > b 



Ucds= , x < b 



Adsorbing surface 



Figure 1 (V. Barsegov and D. Thirumalai) 




Figure 2 (V. Barsegov and D. Thirumalai) 
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Figure 4 (V. Barsegov and D. Thirumalai) 
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Figure 6 (V. Barsegov and D. Thirumalai) 



